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Resumo 

Neste trabalho pretendemos identificar e ilustrar as potencialidades e 
fraquezas do sistema de computagao algebrica Maple, na area do Calculo 
das Variagoes: uma area classica da Matematica que estuda os metodos 
que permitem encontrar valores maximos e mmimos de funcionais do 
tipo integral. Os problemas variacionais sao normalmente resolvidos por 
recurso as condicoes necessarias de Euler-Lagrange, que sao, em geral, 
equagoes diferenciais de segunda ordem, nao lineares e de diffcil resolugao. 
Neste trabalho mostramos como o sistema de computagao algebrica Maple 
pode ser util na determinagao e resolugao das referidas equagoes. Vamos 
tambem apresentar a resolugao do problema classico da braquistocrona, 
sob o ponto de vista do Calculo das VariagSes e do Maple, e uma refor- 
mulagao do problema restringindo a classe de fungSes admissfveis. 
Palavras-chave: sistemas de computagao algebrica, Maple, calculo das 
variagoes, equagoes de Euler-Lagrange, problema da braquistocrona. 

Abstract 

It is the aim of this work to identify and illustrate the potential and 
weaknesses of the computer algebra system Maple in the area of the Calcu- 
lus of Variations: a classical area of mathematics that studies the methods 
for finding maximum and minimum values of functionals of integral type. 
Variational problems are usually solved with the help of the necessary 
optimality conditions of Euler-Lagrange, which are, in general, nonlinear 
and difficult second order differential equations to be solved. We show 
how the computer algebra system Maple can be useful in the determina- 
tion and resolution of these equations. We will also present the solution 
to the celebrated brachistochrone problem from the point of view of the 
calculus of variations and the Maple system, and a reformulation of the 
classical problem obtained by restricting the class of admissible functions. 
Key-words: computer algebra systems, Maple, calculus of variations, 
Euler-Lagrange equations, brachistochrone problem. 
2000 Mathematics Subject Classification: 49-01, 49-04, 49K05. 
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1 Introducao 



O Maple faz parte de lima fanrilia de ambientes computacionais apelidados de 
Sistemas de Computagao Algebrica. A computagao algebrica, tambem chamada 
de computagao simbolica, e uma area de investigacao moderna, que surgiu 
na segunda metade do seculo XX. O Maple e uma ferramenta matematica 
muito poderosa que permite realizar uma miriade de calculos simbolicos. In- 
clui um enorme mimero de comandos, disponiveis em varios packages, os quais 
nos permitem trabalhar em areas como: Calculo das Variagoes, Algebra Lin- 
ear, Equagoes Diferenciais, Geometria, Logica, Estatistica, entre outras. Inclui 
tambem a sua propria linguagem de programagao de alto m'vel que permite, ao 
utilizador, definir os seus proprios comandos. 

O Calculo das Variagoes e uma area classica da matematica, com mais de tres 
seculos de idade e que continua extremamente activa em pleno seculo vinte e um. 
Na sua essentia o objectivo do calculo das variagoes e encontrar um caminho, 
uma curva ou uma supcrfi'cie, para os quais uma determinada funcional tern um 
valor nn'nimo ou maximo. 

Nota tecnica. Neste trabalho usamos a versao 9.5 do Maple. Ao usar uma 
versao anterior ou posterior a 9.5, o leitor podera obter os resultados numa forma 
nao completamcnte identica a aqui apresentada. No essential tudo se mantera, 
no entanto, igual (verificamos isso mesmo usando, a posteriori, a versao 11 do 
Maple). A data presente a versao mais recente do Maple e a versao 12 (cf. 
: //www.maplesof t . com| ). 



2 Calculo das Variagoes 

O problema fundamental do calculo das variagoes consiste na determinagao de 
um extremante - um minimizante ou um maximizante - para uma funcional 
integral que depende da escolha de uma fungao pertencente a uma determinada 
classe de fungoes, em particular fungoes cujos valores nos extremos de um dado 
intervalo real fechado sao fixos. 

Sejam a, b, a < b, A, B numeros reais e L : [a,i)]xKxl — ► M. uma fungao 
continua e diferenciavel com continuidade. Considera-se o problema seguinte: 

.%(•)]=/ L(x,y(x),y'(x))dx — > min, y(a) = A, y(b) = B, (1) 

J a 

onde o minimizante se procura no conjunto das fungoes duas vezes continua- 
mente diferenciaveis no intervalo [a, b}: y(-) £ C 2 ([a, b]). A este problema da-se 
o nome de Problema Fundamental do Calculo das Variagoes. 

Os problemas deste tipo sao normalmente resolvidos recorrendo a equagao 
diferencial de Euler-Lagrange: 

^(x, y(x), y\x)) - £-^(x, y(x), y\x)) = . (2) 
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Definigao 1. As solugoes das equagoes de Euler-Lagrange da-se o nome de 
extremais. 

De um modo geral, a equagao diferencial de Euler-Lagrange e uma equagao 
nao linear, de segunda ordem (ou de ordem superior, quando os problcmas varia- 
cionais envolvem derivadas de ordem superior a um) de diffcil resolucao. 

O package VariationalCalculus do Maple proporciona um conjunto de 
comandos para a resolugao de problemas do Calculo das Variagoes. Para utilizar 
este package temos de o "carregar" , executando para isso o seguinte comando: 

> with (VariationalCalculus) : 

Neste package temos dispom'vel a fungao, 

EulerLagrange (L , x , y (x) ) 

que devolve um conjunto com as equagoes de Euler-Lagrange para uma funcional 
do tipo ([T]). A expressao algebrica L representa o Lagrangiano, x a variavel 
independente e y{x) uma expressao algebrica desconhecida. Esta fungao do 
Maple permite lidar com funcionais com mais de uma variavel dependente: nesse 
caso y(x) £ 1" e representado, em Maple, por uma lista. 

Vejamos uma aplicagao desta fungao Maple no exemplo que se segue: 

Exemplo 1. Encontrar a extremal de Euler-Lagrange y(-) associada a funcional 

J[y{.)}= f{l2x y {x)~(y\x)f)dx 
quando sujeita as condigoes de fronteira 

y(-l) = l e 2/(0) = 0. (3) 
Comegamos por carregar o package em memoria e por defmir o Lagrangiano 

L. 

> with(VariationalCalculus) : 

> L : = (12*x*y(x)-diff (y(x) ,x)~2) ; 



L := 12 xy (x) — 




As respectivas equagoes de Euler-Lagrange sao facilmente obtidas por intermedio 
do Maple. 

Utilizando o comando EulerLagrange obtemos um conjunto que contem a 
equagao de Euler-Lagrange: 

> eqEL := EulerLagrange (L,x,y(x) ) ; 
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eqEL := ^12x + 2—y(x)j 

De notar que, por defeito, quando o segundo membro da equagao e 0, o Maple 
devolve a equagao sob a forma de expressao. 

As extremais sao obtidas resolvendo esta equagao diferencial com as condigocs 
iniciais ([3]) dadas. Recorremos ao comando Maple dsolve: 

> dsolve ({op (eqEL) ,y(-l)=l,y(0)=0>, y(x)); 

y (x) = -x 3 

O comando op foi utilizado para extrair o elemento do conjunto eqEL. 

Em alguns casos particulares a equagao de Euler-Lagrange admite primeiros 
integrals, isto e, quantidades $>(x, y(x), y'(x)) que sao constantes ao longo de to- 
das as solugoes y(x), x G [a, b], da equagao de Euler-Lagrange, permitindo, em 
princfpio, reduzir a ordem da mesma e, consequentemente, simplificar o prob- 
lema. 



O Maple devolve estes primeiros integrals em conjunto com a equagao de 
Euler-Lagrange sempre que e capaz de os encontrar. 

Em |4J consideramos cinco casos de primeiros integrals. Aqui veremos dois 
desses casos: o integral de momenta (no exemplo que se segue) e o integral de 
energia (aquando a resolugao do Problema da Braquistacrona). 

Consideremos o seguinte exemplo: determinar as extremais da funcional 

J[y{x)} = J y'{x){l + x 2 y'(x))dx 

sujeita as condigoes y(l) = 3 e y{2) = 5. 

Como L = y'(x)(l + x 2 y'(x)) nao depende explicitamente de y(x), entao a 
equagao de Euler-Lagrange tern a forma 

e(^W))-o w 

pelo que admite o integral de momento 

dL 

•F7 ( x iV'( x )) = const. (5) 
oy' 

O comando Maple EulerLagrange devolve neste caso urn conjunto com dois 
elementos: a equagao de Euler-Lagrange propriamente dita e o integral de mo- 
mento. Tendo em conta que num conjunto nao importa a ordem de ocorrencia 
dos elementos, nao devemos assumir que os conjuntos devolvidos pelo Maple 
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obedegam ou mantcnham determinada ordem. De maneira a idcntificarmos 
de imediato a ordem dos elementos devolvidos podemos, atraves do comando 
convert, obter uma lista - uma sequencia ordenada de elementos entre parentesis 
rectos: 

> eqEL : = convert ( 

EulerLagrange(diff (y(x) ,x)*(l+x~2*dif f (y(x) ,x)) ,x,y(x)) , 
list) ; 



eqEL := 



~ Ax ic y {x) 2 x2 j^ v {x) ' 1 + 2x2 i y {x) = Ki 



O Maple devolveu uma lista. O primeiro elemento da lista e a equagao de 
Euler-Lagrange abaixo designada por eqEL [1] ; o segundo elemento e o integral 
de momento, abaixo designado por eqEL [2] . Aplicando o comando dsolve a 
equagao de Euler-Lagrange encontramos facilmente a extremal: 

> dsolve ({eqEL [1] ,y (1) =3 , y (2) =5} ,y (x) ) ; 

y(x) = 7-- 

x 

Aplicando o comando dsolve ao integral de momento, nada e devolvido. 

> dsolve ({eqEL [2] ,y(l)=3,y(2)=5},y(x) ) ; 

O Maple fica provavelmente "baralhado" a avaliar e a determinar as constantes. 
Em vez de introduzir as condicoes iniciais no comando dsolve, vamos avalia-las 
mais tarde. 

> sol := dsolve (eqEL [2] ,y(x)) ; 

sol := y (x) = _li±A^± + _ci 

x 

Substituindo os valores de x e y(x), respeitantes as condigoes iniciais, na 
solugao geral obtida, sol, obtemos duas equacoes com duas incognitas. 

> soli := subs({y(x)=3,x=l>,sol) ; 



soh :=3=^-^K 1 +_Cl 

> sol2 := subs({y(x)=5,x=2>,sol) ; 

sol 2 :=5 = \-\k 1 +_C1 
4 4 

Rcsolvemos o sistema de equagoes por intermedio do comando solve: 

> const := solve({soll,sol2},{K[l] ,_C1» ; 

const := {_C1 =7,K 1 = 9} 

Finalmcntc, substitufmos o valor das constantes na solugao geral, sol, e obtemos 
a extremal pretendida: 
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> subs ({op (const)}, sol) ; 

y(x) = 7 

x 

Notamos que se derivarmos em ordem a x o integral de momento obtemos a 
equagao de Euler-Lagrange. Para vermos isto com o Maple utilizamos o comando 
diff: 

> diff (eqEL[2] ,x) ; 

ix^-y(x) + 2x 2 ^y(x) = 

Podemos dizer que do ponto de vista matematico a dificuldade de determi- 
nar as extremais por intermedio da equagao de Euler-Lagrange ou do primeiro 
integral e a mesma, mas do ponto de vista do Maple nao: embora equivalentes, 
o Maple nao resolve da mesma forma a equagao de Euler-Lagrange e o primeiro 
integral. Embora o processo utilizando o integral envolva mais calculos, obtive- 
mos, como esperado, a mesma solugao. 



Para problemas deste tipo, o calculo das extremais, sob o ponto de vista do 
Maple, revelou-se mais simples recorrendo a equagao de Euler-Lagrange do que 
utilizando apenas os respectivos primeiros integrals. No entanto so encontramos 
uma solugao satisfatoria para o problema da Braquistocrona atraves da resolugao 
a partir do primeiro integral. O leitor interessado na tematica dos primeiros 
integrals no calculo das variagoes podera ver [2] . 



3 O Problema da Braquistocrona 

Um dos primeirosQ e, sem diivida, o mais famoso problema associado ao desen- 
volvimento da teoria matematica do calculo das variagoes e, como ja foi referido, 
o problema da Braquistocrona (para uma descrigao historica e extremamente 
interessante deste problema recomcndamos o excelente artigo 

Dados dots pontos A e B num piano vertical, qual o caminho APB que a 
particula raovel P atravessa em tempo minimo, assumindo que a sua aceleragao 
e apenas devida a gravidade? 

O problema da Braquistocrona e formulado matematicamente como se segue: 



T[y(-)} = 


i 


L i 








2/o, 





1 + (y'(x)) 2 J . 

—^—dx — > mm (6) 

Vo - y(x) 

= 2/1, y € C 2 (xq, xi) 



1 outro problema famoso que esteve na origem do calculo das variagoes e do controlo 
optimo e o problema de Newton da resistencia minima, colocado e resolvido por Isaac Newton 
em 1686 8 . Para desenvolvimentos recentes deste problema veja-se [6] [9]. 
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sendo T o tempo necessario para a particula dcslizar da posigao inicial A(xq, j/o) 
ate a posigao final B(xi, y\). 

A solugao deste problema e um arco de uma cicloide gerada pelo movimento 
de um ponto fixo localizado numa circunferencia, de diametro a, que rola sobre 
a recta y — yg e cujas equagoes parametricas sao dadas por 

' x = x + § (0- sin(0)) 

(7) 

,y = W)-f(i-cos(fl)) 

com 6*o < 9 < 9i, onde 6q e 9\ sao os valores de 9 nos pontos A(xo,yo) e 
B{x u yi). 

O par de equagoes 

a (9i — sirx(^i)) = 2(xx — xq) 

(8) 

a(l- cos(6»i)) = — 2(yi - y ) 

permite determinar valores unicos para as constantes a e 9\ (0 < 9\ < 2ir) em 
fungao dos valores dados xq, Xt, yo e y±. 

3.1 Solugao da equagao de Euler-Lagrange no caso da Braquistocrona 

No Problema da Braquistocrona © , ignorando o factor constante -^= que nao 
altera em nada o minimizante, a fungao ?/(•) que minimiza T sera encontrada 
atraves da equagao de Euler-Lagrange ([2]) com 



Como L nao depende explicitamente de x, a equagao de Euler-Lagrange admitc 
o integral de energia 

dL 

L{y{x),y'{x)) - y'(x) — (y(x),y'(x)) = const, 

ou seja, 

dL 

y'( x )-Q^j(y( x )> y'( x )) - L (y( x ), y'i x )) = const. 

Tendo por base [T] e [TU], segue-se uma analise do problema recorrendo 
aqui ao sistema de computagao algebrica Maple. 

> L := (1 + diff (y(x) ,x)~2)~(l/2)/(y0-y(x))~(l/2) ; 



\Jyo - y (x) 
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> eqEL := convert (simplif y(EulerLagrange(L,x,y(x) )) ,list) ; 



eqEL : 



+ W) 2 ~ 2 (*,y (xj) yo + 2 (£y (x)) y (x) 
(w-!/rf(l + (|»W)f 
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#1 



Neste caso, a equagao de Euler-Lagrange possui uma certa complexidade e 
o Maple nao consegue obter uma solugao satisfatoria. 

Nao sendo vantajoso trabalhar com a equagao de Euler-Lagrange, vamos 
trabalhar com o integral de energia e determinar com o Maple uma forma 
parametrica das extremais para o problema da Braquistocrona. Ao aplicar o co- 
mando dsolve ao integral de energia, com a opgao parametric, o Maple devolve 
as extremais para o problema na forma parametrica: 

> sol := dsolve (eqEL [2] ,y(x) .parametric) ; 



As expressoes matematicas de x e y sao dadas em fungao do parametro _T. 
A constante -CI e a segunda constante de integragao. 

A representagao encontrada pelo Maple e menos elegante do que a repre- 
scntagao parametrica usual da cicloide, atras descrita. Conseguimos no entanto 
obter ([7]), via Maple, atraves de alguns artiffcios. 

Uma nova fungao 9 = 6{x) e introduzida atraves da seguinte relagao: 
> _T := -cot(theta/2) ; 




yo K 1 2 + y K^.T 2 - 1 




z(_T) 



_T + arctan (_T) + arctan (_T) _T 2 + _C1 K x 2 + .CI K^.T 2 
K x 2 (1 + _T 2 ) 



(10) 



T 




> K[l] := l/sqrt(a) ; 
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* - 4. 

O segundo membro (right-hand side) da segunda equagao devolvida em sol 
e dado, em Maple, por rhs(sol [2] )): 

> x := rhs(sol[2]) ; 



x := 



(- cot (|) - f + arccof (cot (|)) + (-§ + arccoi (cot (f ))) (cot (f )) 2 ) a 



l + (cot(f)) 2 



_Ci+_Ci (cot(f)) 2 
l + (cot(|)) 2 

A expressao obtida pelo Maple e claramente possfvel de ser simplificada, bas- 
tando, para isso, substituir arccot (cot (|)) por |: 

> x := simplify(subs(arccot(cot(theta/2))=theta/2,x)) ; 

aj:= _£ sin ( ( 9) + £0_Z[ a + _ C rj 

Recordando que 6* = quando x — x 0} temos _Ci = x + §a: 

> x := subs(_Cl=xO+(Pi/2)*a,x) ; 

a; := -7; sin (8) + ^8 + xO 

Seleccionando o segundo membro da primeira equagao devolvida em sol (rep- 
resentado, em Maple, por rhs (sol [1] ) ) e aplicando o comando simplify, obte- 
mos a forma final desejada para y: 

> y := simplify (rhs (sol[l] )) ; 

a a 
y:=- cos (6) - - + y 

Com a forma parametrica de x e y determinada, o tempo de descida mfnimo, 
T, pode ser facilmente calculado em Maple. Derivando x e y em ordem a 9, 
obtemos: 

> x'(theta) := dif f (x,theta) ; 

x'{6) :=£(l-cos(0)) 

> y'(theta) := diff (y.theta) ; 
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1/(0) :=--sin(0) 

De forma a calcular 

+ dx{6) 

fazemos: 

> f := simplify(sqrt(l+(y'(theta)/x'(theta))"2)*x'(theta)/sqrt(y0-y)) ; 

/ := - (-1 + cos (0)) ^-(-l + cosWr 1 a === 

V -a (-1 + cos (0)) 

Vamos simplificar a expressao / assumindo 9 > 0. Designaremos a expressao 
devolvida por "integranda" . 

> integranda := simplify (f) assuming theta > 0; 

integranda := \fa 

O tempo mmimo de descida, T, calcula-se integrando a "integranda" de 
9 = ate 6>i: 

> Tthetal := (l/sqrt(2*g))*int(integranda,theta=0. .thetal) ; 

Tthetal := ^1 

Os valores das constantes a e 9i dependem das coordenadas dos pontos A e B. 
Para a constante gravitacional g tomemos o valor aproximado 9.8: 

> g := 9.8: 

3.2 Exemplos 

Exemplo 2. Calculemos o tempo de descida T para o caso em que A(0, 2) e 
£>(3, 1) c traccmos a curva de descida que conduz a esse tempo T. 

Tal como referido anteriormente, o par de equagoes ijHJl permite determinar 
valores unicos para as constantes a e 9\ (0 < 9\ < 2ir) em fungao das coor- 
denadas xq, Do e xi, y\ dos pontos A e B, respectivamente. Vamos definir 
fungoes apropriadas em Maple que nos permitem resolver o problema de forma 
generica. Ilustramos as nossas fungoes Maple para o caso particular em que 
(xo,j/o) = (0,2) e {xi,yi) = (3,1). Determinemos entao os valores de a e 9\ 
para A(0,2) e 5(3,1): 

> eql := (xO.xl) -> a*(thetal-sin(thetal))=2*(xl-x0) ; 

eql := (af0,a;l) -> a(91 - sin(01)) = 2x1 - 2x0 

> eq2 := (y0,yl) -> a*(l-cos(thetal))=-2*(yl-y0) ; 
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eq2 := {yO, yl) -» o(l - cos(01)) = -2t/l + 2y0 

> sol := (xO,yO,xl,yl) -> f solve({eql(xO,xl) ,eq2(y0,yl)},{a,thetal}) ; 

sol := (xO, yO, xl, yl) — ► fsolve(eql(xQ, xl), eq2(y0, yl), a, 91) 

Os valores de a e #1 para o nosso exemplo sao obtidos chamando a fungao sol 
com os valores de xq = 0, yo = 2, xi = 3, y\ = 1: 

> sol(0,2,3,l) ; 

{a = 1.239374053 , 91 = 4.051628024} 

Para obter os valores de a e 0\ separadamentc, defmimos as funcocs va ( "valor 
de a") e vthetal ("valor de &i"): 

> va := (x0,y0,xl,yl) -> eval(a,sol(xO,yO,xl,yl)) : 

> va(0,2,3,l); 

1.239374053 

> vthetal := (x0,y0,xl,yl) -> eval(thetal,sol(xO,yO,xl,yl)) : 

> vthetal (0,2, 3,1) ; 

4.051628024 

Dcfinimos agora em Maple a funcao Tminimo(xo, yo, %i,yi) que nos permite 
calcular o tempo de descida mi'nimo T c cxprimimos a extremal parametrica- 
mente atraves de soIx(xq, yo, x\, y\) e soly(xo 1 yo,xi 1 yi): 

> Tminimo := (x0,y0,xl,yl) -> evalf (eval (Tthetal , sol(x0 ,y0 ,xl ,yl) ) ) : 

> T : = Tminimo (0,2, 3, 1) ; 

T := 1.018832361 

> solx := (xiO,yO,xl ,yl) -> eval (x, sol (xi0,y0,xl,yl) union {x0 = xiO}) : 

> solx(0,2,3,l) ; 

0.6196870265 6 - 0.6196870265 sin (6) 

> soly := (x0,yi0,xl ,yl) -> subs(sol(xO,yiO,xl,yl) union {yO = yiO},y): 

> soly(0,2,3,l) ; 

1.380312974 + 0.6196870265 cos (6) 
Ao longo da curva defmida parametricamente por 

x = 0.6196870265 8 - 0.6196870265 sin (9) , 
y = 1.380312974 + 0.6196870265 cos (9) , 

a partfcula demora aproximadamente 1.02 scgundos a deslizar de A ate B. Este 
6 o tempo mmimo aproximado entre estes dois pontos. O leitor e convidado a 
usar as funcocs Maple que acabamos de definir nos seus proprios exemplos. 
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Recorrendo ao Maple, calculamos tambem facilmente o tempo de descida da 
partfcula ao longo de outras curvas que unam A a. B. 

> L := (y,yO) -> sqrt ( (1+dif f (y (xi) ,xi) ~2) I (yO-y (xi) ) ) ; 



1 



\ yO -y (?) 



> Tempo := (y,xO,xl,yO) -> evalf ((l/sqrt(2*g))*int(L(y,yO) ,xi=xO. .xl)) ; 



Tempo := (y, xO, xl, yO) — > evalf —7= / L(y, yO)dt; 



xO 

Por exemplo, se determinarmos o tempo de descida, Ti, ao longo da recta 
2/i definida por A e B, verificamos que e superior a T. 

> yl := x -> -(l/3)*x+2; 

yl := x — > -- + 2 
y 3 -r 

> Tl := Tempo (y 1,0, 3, 2) ; 

Tl := 1.428571428 

Seguindo a proposta feita em |5J (embora com outros valores), podemos 
tambem determinar o tempo de descida, T2, ao longo do arco da circunferencia 
com centro (3,6) e que passa pelos pontos A(0, 2) e 5(3, 1) e verificar que e maior 
do que T. Esta circunferencia tern por equagao reduzida: (x-3) 2 + (y-6) 2 = 25. 

Consideremos a semicircunferencia definida por yi'. 

> y2 := x -> 6-(16-x~2+6*x)~(l/2); 



y2 := x — > 6 — \/ 16 — x 2 + 6 x 

Calculamos facilmente o correspondente tempo Ti por intermedio da nossa 
funcao Maple: 

> T2 := Tempo (y2, 0,3, 2) ; 

T2 := 1.151743820 

O Maple permite criar representacpes graficas de varias curvas, separada- 
mente ou em simultaneo. Atraves do comando display, dispomvel no package 
plots do Maple, conseguimos obter uma representacao grafica das tres curvas 
acima consideradas: cicloide, recta e arco de circunferencia. Para podermos uti- 
lizar este comando temos de carregar o package no qual ele se encontra inserido: 

> with(plots) : 

Defimndo pi, P2 e P3 como sendo as "estruturas plot" 
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> pi := plot( [solx(0, 2,3,1) ,soly(0, 2,3,1) , 

theta=0. .vthetal(0,2,3,l)] , 
scaling=constrained , thickness=2 , 
labels= ["x" , "y"] , color=blue) : 

> p2 := plot(yKxl) ,xl=0. .3, thickness=2, labels=["x" , "y"] , 

scaling=constrained, color=black) : 

> p3 := plot(y2(x2) ,x2=0. .3, thickness=2, labels=["x" , "y"] , 

scaling=constrained, color=red) : 

obtemos entao uma representacao grafica das tres curvas acima consideradas 
(cicloide, recta e arco da circunferencia) atraves do seguinte comando display: 

> display(pl,p2,p3) ; 




Embora a cicloide seja a curva de maior comprimento, entre as tres comparadas, 
e a curva do tempo mi'nimo. 

Exemplo 3. Calculemos o tempo de descida minimo T para o caso em que A(l, 3) 
e -B(15, 1) e tracemos a curva. 

Para este problema temos 6\ dado por 

> vthetaK 1,3, 15,1) ; 

4.875635855 

enquanto a toma o valor de 

> va(l,3,15,l) ; 
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4.776249228 

O tempo de descida mmimo e inferior a 2.41 segundos 

> TminimoC 1,3, 15,1) ; 

2.406837209 

sendo alcangado atraves da curva descrita parametricamente por 

> x := solx(l,3,15,l) ; 

x := 1 + 2.388124614 9 - 2.388124614 sin (9) 

> y : = soly(l,3,15,l) ; 

y := 0.611875386 + 2.388124614 cos (6) 

Podemos tracar facilmente esta curva em Maple usando o comando plot 
e fazendo variar o 9 de ate 9\. Matematicamente, a curva e parte de uma 
cicloide invertida. 

> plot( [x,y,theta=0. . vthetal (1 ,3 , 15 , 1)] , 

scaling=constrained, thickness=2, labels= ["x" , "y"] ) ; 




4 Uma Variante do Problema da Braquistocrona 

Para alem do problema publicado por John Bernoulli apresentado em muitas 
variantes do problema da Braquistocrona tern sido colocadas e resolvidas ao 
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longo dos tempos: Braquistocrona atraves do planeta Terra 10J ; problema sob 
a presenca de atrito j3]; etc. 

Na solucao classica foi provado que a Braquistocrona une quaisquer dois 
pontos do piano vertical. Isto porque se assumiu a Terra como inexistente: 
apenas o campo gravitacional existia. 

De acordo com [7], na reformulagao do problema classico, que se segue, va- 
mos impor que a linha da Braquistocrona nao passe para valores negativos de 
Oy. Isto implica, obviamente (cf., por exemplo, a figura no final de £)3.2p , que a 
solugao geral do problema sera diferente. 

Sejam A(0,1) e B(b,0), b > 0. Ignorando o factor constante recon- 
sideramos o problema da Braquistocrona restringindo a classe das fungoes ad- 
missfveis: 

T[y(-)} = i/ ^j^ dx — > min, yeS, (11) 

com 

S={y(x):y(0) = l, 2/(1) = 6, y" > 0, < y(x) < y (x) , yeC 2 (0,6)} 
onde yo(x) e a recta que passa pelos pontos A e B e tern por equacao 

yo(x) = 1 - -. 

Nesta reformulagao do problema nem todos os pontos podem ser unidos por 
uma extremal da funcional (fTTj) como se prova em [J] utilizado um novo metodo 
de integracao da equagao de Euler-Lagrange. 

Em [4], e explorada tambem a seguinte questao: 

Sera verdade que 

T(y)<To:=T(y Q ) = 2Vl + P, Vy E S ? 

Servindo-nos das potencialidades do Maple (em particular, recorremos a 
fungao Thiele Interpolation dispomvel no package CurveFitting), encontramos 
uma fungao y € S que verifica a desigualdade 

T(y) > T . 

O leitor interessado encontrara os detalhes em [4] . 
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